Weak localization correction to the FS interface resistance 
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The classical resistance of a contact between a ferromagnet (F) and a superconductor (S) 
acquires an additional contact term as compared to the contact between a ferromagnet and a normal 
metal. The necessity to match spin-polarized current in a ferromagnet to spin-less current in the 
superconductor results in the accumulation of non-equilibrium polarization near the F/S interface. In 
the present work, we show that the weak localization correction to the classical diffusion coefficient, 
5D, is dependent on the degree of polarization, with majority spins more likely to be reflected from 
the interface than minority spins. Taking into account the change in the spin polarized particle 
distribution in the F wire arising from 5D, we calculate the weak localization correction to the F/S 
contact resistance. 



PACS numbers; 73.23.-b, 72.10.Bg, 74.80.Dm, 



I. INTRODUCTION 

Novel phenomena in mesoscopic systems have been the 
subject of inteBse theoretical and experimental interest 
for many yearsoH. Observed effects such as weak local- 
ization and universal conductance fluctuations are due 
to the quantum interference of electrons at low tempera- 
tures. An important branch of mesoscopic physics is that 
of hybrid nanostructures where the influence of a super- 
conductoij-||(S) on a phase coherent normal (N) region 
is studiecila. At low temperatures, the role of Andreev 
reflectionEl is important, whereby particles in the N region 
with excitation energies e lower than the superconduct- 
ing gap energy A are reflected from the N/S interface as 
holes. Studies of subgap transport have shown that su- 
perconducting condensate penetration into the N region, 
known as the proximity effect, may substantially change 
the resistance of an N/S circuit. 

Lately, improvements in the microprocessing of metals 
has led to the fabrication of nanostmctures combining 
ferromagnetic (F) and S materialsO^O. The presence of 
a large exchange field, Scx ^ A, suppresses electron-hole 
correlations in ferromagnets so that the role of the prox- 
imity effect is reduced-. However a separate mechanism 
has been predicteclE3~tJ to produce a resistance increase 
in a circuit consisting of a mono-domain F wire connected 
to an S electrode instead of an N electrode. This is a ro- 
bust, classical effect which originates from the need to 
match a spin-polarized current in the ferromagnet to a 
spinless current in the superconductor. 

The present work analyses the weak localization cor- 
rection to the contact F/S resistance. In a semiclassi- 
cal language, weak localization arises from an enhanced 
backscattering caused by the xuiaiitum interference of 
pairs of coherent quasi-particlesEjIlj. Constructive inter- 
ference of a quasi-particle paired with its time reversed 
partner is destroyed by a magnetic field so we assume 
that the width of the F wire, L±, is small, L± < lOOnm, 
in order to limit the influence of the intrinsic magnetic 



field. Weak localization is affected by the boundary con- 
ditions of the system. Particles that escape into an N 
electrode suffer dephasing which reduces the return prob- 
ability. However, multiple Andreev reflections from an S 
electrode may enable a particle to return coherently. For 
a polydomain F wire, with no net average polarization, 
the return probability [Of an F/S system is greater than 
that of an F/N systemtS. 

In general, the spin channels have different conductiv- 
ities, aa = e^VaDa, whcrc i/a and Da are the density of 
states and the classical diffusion coefficient for electrons 
in the spin state a, a — {'], I)- The degree of spin polar- 
ization, ^, is defined as C = ((T| ~ cx) / (cj -I- (T|) and we 
calculate the return probability in an F/S system with 
an arbitrary degree of polarization, < C < 1. In ad- 
dition we introduce a finite spin relaxation length, L^, 
shorter than the length, L, of the F wire, Lg ^ L. As 
the spin polarization increases, the majority spin chan- 
nel experiences a further increase in return probability 
(as compared to the F/N case), whereas the minority 
spin channel suffers a reduction. Majority carriers at 
the superconducting interface cannot find minority car- 
rier states to Andreev reflect into and are more likely 
to be reflected normally whereas minority carriers at the 
interface have an enhanced probability to escape. For 
C ^ 1 the reservoir appears to be totally insulating to 
majority spins but normal as far as minority spins are 
concerned. 

After calculating the correction to the classical diffu- 
sion coefficient that is due to the enhanced return prob- 
ability, we find the corresponding change in the spin po- 
larized particle distribution and we determine the weak 
localization correction to the classical resistance. Gener- 
ally this is a sum of a bulk term that is independent of the 



state of the reservoir and a contact term, Sr, 



FS{FN) 



that 



depends on whether the reservoir is in the superconduct- 
ing or the normal state. We find that Sr^^ > 6r^^ , so 
that a change in the state of the reservoir from N to S re- 
sults in an increase of the weak localization contribution 
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to the resistance, 



6r: 
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Sr, 



FN 



(1) 



where Ra = 1/((T|+o'|)L'[^^ is the resistance of a cube of 
length equal to the width of the wire L±. The combina- 
tion RaLs/L± corresponds to the resistance of a piece of 
ferromagnet of length Lg. We evaluate the prefactor, 7, 
numerically and find that it is almost constant, 7 « 1/2, 
for all experimentally relevant values of spin polarization. 



II. CLASSICAL RESISTANCE 

We begin by briefly describing the approach of 
Refs. ^0 which enables one to calculate the classical 
contact resistance of an F/S interface and we refer the 
reader to those papers for further details. We consider a 
single domain, disordered ferromagnetic wire of length L 
as depicted in Fig. 1. The coordinate x is used to describe 
the position along the wire and throughout the paper we 
consider there to be a normal reservoir at the left hand 
side, X = — L, whereas the reservoir at the right hand 
side, a; = 0, is either N or S. The resistance of the disor- 
dered F wire is determined by solving diffusion equations 
for the isotropic part of the electron distribution function 
na{s, x) = J dflpna{p, x). Using electron-hole symmetry, 
it is possible to consider only the symmetrised function 
Na{e,x) = [na{£,x) + na{—£, x)]/2, where e is the en- 
ergy determined with respect to the chemical potential 
in the S (N) reservoir. The current density in any given 
spin channel is related to the spatial derivative of the 
distribution function as 



Jc 



dedxNa{e,x). 



(2) 



The distribution function obeys a differential equation 
describing diffusion in the ferromagnetic wire, 

DadlNa{e,x) ^ w^ii^[Na{e,x) -Na{e,x)], (3) 

where a = (i,t) for a = (tii)- This equation may be 
expressed conveniently as 



dl[a^N^+a^Ni] = 0, 
{dl - L-^) [TVt - 7VJ = 0. 



(4) 
(5) 



Spin relaxation, which may arise from spin-orbit scatter- 
ing or spin-flip scattering at defects, is described by the 
right-hand side of Eq. (H). One may define an effective 
spin relaxation rate, Lg, as Lj^ = w-^ ^[v'^ / D ^ + v^/D^] 
which accounts for the relaxation of the difference be- 
tween the spin channels. In a similar way, the spin re- 
laxation length in an individual channel, La, is defined 



as L^"^ = w^[U^/Da such that L^"^ 
equivalent expression is = Ll{aa 
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The above equations should be complemented by four 
boundary conditions, two at each end of the ferromag- 
netic wire. The boundary conditions at the left hand 
side of the wire are given by the equilibrium distribution 
of electrons in the normal reservoir, 

iV„(e, -i) = i [nris - eV) + nri-e - eV)] , (6) 

such that the distribution of up and down spins is equal. 
The boundary conditions at the right hand side of the 
ferromagnetic wire at the superconducting reservoir have 
been discussed in detail in Refs. For quasiparticles 

with energies below the superconducting gap they may 
be written as 

a^d.N^\,= a^d.N^\^, (7) 
N^{£, 0) -f 7Vx(£, 0) = [7iT{e) + nri-e)] . (8) 

The first of these, Eq. (^), describes Andreev reflection 
such that the spin-up current must equal the spin-down 
current and the second, Eq. (^, states that the sum of 
the distributions is equal to the equilibrium value. On 
solving the above differential equations with the appro- 
priate boundary conditions it is possible to calculate the 
current density in the wire. In particular, the spatial 
derivative of the distribution is 



dxNa{,e,x) 



(1 - 2iVi) 



2rL 



(Ja) cosh[(a; + L)/Ls 



where 



r = 1 



2a a 



cosh(L/Ls) 



tanh{L/Ls), 



(9) 



(10) 



and Nl = Na{e,-L). 

The effect of the superconducting reservoir, for strong 
spin relaxation L ^ L^, is expressed in terms of a con- 
tact resistance r^^ so that the total classical resistance, 
Rci , may be written as a sum of resistances in series 



L 

Rcl — 7 — Rn 
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Rl 



Ls {(7] - criY 



(11) 



(12) 



For a normal reservoir on the right hand side, the bound- 
ary conditions are similar to those in Eq. (^), namely 
Na{£, 0) = [nrie) + nT{—e)\ /2. In this case, there is no 
contact resistance. A fall in temperature, precipitating a 
change in state of the reservoir from N to S, would there- 
fore result in a resistance increase of the circuit. This is 
a robust classical effect which originates from the need 
to match a spin-polarized current in the ferromagnet to 
a spinless current in the superconductor. It is specific 
to mono-domain wires since in poly-domain wires the 
tranport properties of spin-up and spin-down particles 
do not differ and, on the average, the current is not spin- 
polarized. 
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III. WEAK LOCALIZATION CORRECTION TO 
THE CONTACT RESISTANCE 

The weak localization correction to the classical diffu- 
sion coefficient, SDa, is related to a two partidjeJ^reen's 
function known as the Cooperon propagatoill3il3. The 
Cooperon consists of a quasi-particle following a diffu- 
sive path that interferes with a quasi-particle travers- 
ing the same path in the opposite direction. When the 
particles follow a closed path, the interference results in 
an enhanced return probability. In a ferromagnet, sin- 
glet and Sz = triplet Cooperons are suppressed since 
the Fermi momentum for spin-up and spin-down parti- 
cles is different, producing a difference in the phase ac- 
cumulated by the particles along the path. However, 
there is no such suppression of the phase correlation in 
the triplet channel where same spin particles are paired. 
Thus we assume that the return probability is given by 
the triplet CoopercuipM-opagator, Caaix, x'), for equal 
spatial coordinates 



DO 



(13) 



which obeys a diffusion equation in the disordered ferro- 
magnet, 

{-VaDadl + i^cT-'^) Caa{x, x') ^d{x~x'). (14) 

Here the spin relaxation time, Tq, is related to the spin 
relaxation length in the spin channel a, Ta = L'^/Da- 

We stress that the existence of a Cooperon is not a re- 
sult of the proximity effect as we do not consider triplet 
pairing induced in the ferromagnet by the presence of 
the superconductor. Nevertheless, weak localization is 
affected by the boundary conditions of the system. A 
schematic of the change in the boundary conditions of 
the Cooperon at the F/S interface as compared to the 
F/N interface is shown in Fig. 2. A process which pairs 
two spin-up quasi-particles is depicted in Fig. 2a. When 
the reservoir is in the N state, particles escaping into it 
suffer dephasing, thus destroying the Cooperon. At the 
left hand reservoir (x = —L), which is in the N state, we 
therefore have 



CTT(-i,x') = Cu(-i,x') = 0. 



(15) 



When the right hand reservoir {x = 0) is in the S state 
we apply 

CTT(0,.T')=c'^Cu(0,a;'), (16) 

cTt d^C^^{x,x')\^^^ = -aie'^ d^Cnix, x')\^^^ , (17) 

where x is a phase accumulated on reflection from the 
superconductor. These boundary conditions, similar to 
those for the density Eqs.(0,||), account for Andreev re- 
flection. Multiple Andreev reflections at an S electrode 
may enable a particle to return coherently to the orig- 
inal place in the F wire with the same spin polariza- 
tion. Fig. 2b shows a process whereby a pair of spin-up 



quasi-particles are Andreev reflected as spin-down quasi- 
holes which are subsequently Andreev reflected as spin- 
up quasi-particles. The calculation of 6Da by solving 
the Cooperon diffusion equation with the above bound- 
ary conditions is described in Appendix 

The weak localization correction to the diffusion coef- 
ficient, 6 Da, leads to a correction to the current density, 
Sja = e^VaSDa de dxNa{e,x). However, 5Da also 
produces a correction to the particle distribution func- 
tion. We take this into account perturbatively by assum- 
ing that the correction, (5-Dq., to the classical diffusion 
coefficient, I?", is small and that there is a correspond- 
ing small correction, 6Na, to the classical particle dis- 
tribution, N^- Expressions for the total diffusion coeffi- 
cient. Da — D^+5Da, and the total particle distribution, 
Na = Na + SNa , are substituted into the previous differ- 
ential equations describing diffusion in the ferromagnet, 
Eqs. (§,|1), and the boundary conditions, Eqs. (§,@jg), to 
provide new equations for SNa- Hence the differential 
equations for the correction to the distribution function, 
SNa, in terms of the correction to the classical diffusion 
coefficient, SDa, are 



5J^ 



(18) 
(19) 



where 5Ja{x) — e^VaSDaix) dxNa{x). The boundary 
conditions at the ferromagnetic reservoir on the left hand 
side are 



5N^{e,-L) = SNi{e,-^L) = 0. 



(20) 



whereas the boundary conditions at the superconducting 
reservoir on the right hand side are 

a^dxdN^l^ ~ a^dxSN^l^ = SJ^ - SJ^, (21) 
SN^{e,0)+SN^{e,0)^0. (22) 

General solutions to the differential equations, 
Eqs. (|l|,|l9|), are 



[5My)+6J^{y)]dy, 
(23) 



and 



cosh ( — I dy. 



(24) 



where U,V,W, and Y are unknown factors to be deter- 
mined by the four boundary conditions. It is thus pos- 
sible to evaluate 6Na in terms of the correction to the 
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diffusion coefficient 
distribution, iV°. 

Taking into account botli SDa and SNa, the weak lo- 
calization correction to the current density summed over 
both spin channels is 



5Da, and the classical part of the gives 7 1 — cr|/((T| + <Ji) ^ 1. This rough estimation 

of the value of 7 at ^ = 1 appears to be in agreement 
with Fig. 3. However, some notes of caution about the 
results for large polarization need to made. Firstly, the 
graphs are shown as a function of ( for fixed L/Ls — 20, 
whereas the spin relaxation lengths, La, in the individual 
spin channels depend on ^. In particular, L| ^ Lg for 
large Secondly, the error in the numerical evaluation 
of the diffusion coefhcient increases rapidly as C — > 1. 
Finally, values of C > 0.8 are unlikely to be realised in 
real materials and, even if they were, the intrinsic mag- 
netic field in such materials would destroy any quantum 
interference. 



/•OC 

/ de [SJ^ + SJi + a^dJN^ + iidJN^] . (25) 
Jo 



This is determined by the factor U in the general solution 
Eq. (p3|). Using the boundary conditions, Eqs. ( pO| , ^ , p2| ) , 
to evaluate U, the weak localization correction to the cur- 
rent density may be expressed as 



de 

Y 



dx J 

SJ] 



[5J^ + SJ^] + 

cosh [{x + L)/Ls 



cosh [L/Ls 



(26) 



The term in brackets contains spatial dependences 
arising from the cosh term and from SJ^aix) = 
e^VaSDaix) d^N^ix), where d^N^^x) depends on 
cosh[(a; -I- L)/Ls], Eq.(^, and dDa{x) is given in Ap- 
pendix^. Taking these terms into account, it is straight- 
forward to perform the integration with respect to the 
coordinate x. The final step is to evaluate 5Da{x) nu- 
merically as detailed in Appendix For a long wire, 
we express the weak localization correction as a sum of 
resistances. 



5R = 5R^ + Srf'^ 



(27) 



where SR^ is the bulk contribution, 5R^ - 
(8e^/7r)i?pi(L| -I- Li)/L\ and 5r^^ is a contact term. 



IV. CONCLUSION 

The weak localization correction to the classical diffu- 
sion coefficient, 5D, is dependent on polarization, with 
majority spins more likely to be reflected from the F/S 
interface than minority spins. Taking into account the 
change in the spin polarized particle distribution in the 
F wire arising from 5D, we found that the weak local- 
ization correction to the contact resistance is related to 
the square of the resistance of a piece of F wire with 
length equal to the spin relaxation length (see Eq. (||)) 
with a numerical prefactor that is almost constant for all 
experimentally relevant values of spin polarization. 
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Ru^] , (28) 



The parameter 7 is plotted as a function of the degree 
of polarization, C, in Fig. 3. The solid line is 7 for 
L/Ls — 20 whereas the dashed line is the contribution 
of the first two terms in Eq. ( p5[ ) only, neglecting the in- 
fluence of density variations. It is worth noting that the 
analytic result for an F/N system is 7 = and for a polyr 
domain wire connected to an S reservoir it is 7 = 1/2113. 
In these cases, 5N — 0. For moderate values of C,, we 
find that there is no large change in 7, 7 «:! 1/2, and 
the influence of density variations is small. However, for 
very large spin polarization, C, ^ 0.8, 7 increases dramat- 
ically and the role of density variations is vital. Drawing 
on the interpretation of the results for the return prob- 
ability mentioned in Appendix we naively estimate 
the resistance of a system with C = 1 by considering the 
majority (up) spin channel to behave as if connected to 
an insulating barrier and the minority channel to behave 
as if connected to a normal barrier. Such a procedure 



APPENDIX A: CORRECTION TO THE 
DIFFUSION COEFFICIENT 

This appendix describes the evaluation of the weak 
localization correction to the classical diffusion coeffi- 
cient, 5Da, Eq. (|T^). In order to calculate the Cooperon, 
Eq. (|l^), we consider spinor eigenvectors with compo- 
nents ^Jna{x) that obey the following diffusion equations: 



{^-v^Didl + v^T^^^ i^niix) = A„-0„x(a;), 
such that the Cooperon Caa is given by 



CoiOL {,X, X^ 



^ \i^na{x)Y 



Xr, 



(Al) 
(A2) 

(A3) 



The momenta in the spin channels, Qna, are related, 
v^D^Ql^ + i^TTf ^ = ^iDiQl^ + v^T-'^ = A„. At the 
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left hand reservoir {x = —L) the Cooperon is zero since 
a particle that escapes into it suffers dephasing, Eq. (15), 

V'„T(-i) = ^ru(-i) = 0. (A4) 

whereas when the reservoir on the right hand {x — 0) is 
a superconductor, the boundary conditions Eqs. ( |l6|jl7| ) 
give 

V'„T(0) = e^'^V«i(0), (A5) 
cr| ^^1pn^{x\^Q = -cTie'^ 9^V'ni(a:)|^^n . (A6) 

The spinor eigenvectors also obey a normalisation con- 
dition, S°^{\il)n^{x)\'^ + \'ipni{x)\'^)dx = 1. We find the 
general solutions to the diffusion equations and use the 
boun dar y conditions at the left hand lead {x = —L), 
Eq. ([A^), to give 



-ipnaix) = N'na Sm [Qnaix + L)] 



(A7) 



On substituting into the boundary conditions on the right 
hand side {x = 0), we find 

MiT sin [Qn]L] = e'^Afni sin [Q„iL] , (A8) 
a-f Q n-fAfni cos [QniL] = -e'^aiQniAfniCos[QniL]. (A9) 

Eliminating the normalisation constants leaves an equa- 
tion for determining the eigenvalues. 



f^lQnicos [QniL] sin [Qn]L] + 

cr|Q„l cos [Qn-fL] sin [QniL] = . 



(AlO) 



We solve this equation numerically to find the eigenval- 
ues and determine their contribution to the Cooperon, 
Eq. @. 

As an example we present results for the calcula- 
tion of the spatially averaged return probability, SDa = 
SDa(x){dx/L), for arbitrary spin polarization in the 
ferromagnet. For long wires, L 3> Ls, we express it as 

La 



L 



(All) 



where SD^ is the bulk contribution, SD^ = 
— SLq/tti^qL^. The numerical coefficient r]a depends on 
the state of the reservoir and the degree of polarization 
in the ferromagnet. We jCompare with analytic results 
obtained in certain limitst^. When the reservoir on the 
right hand side of the wire is in the normal state, the 
boundary conditions are the same as those in the left 
hand reservoir and r/^ = f . For a poly-domain ferromag- 
netic wire connected to an S reservoir, where the classical 
contact resistance Eq. (f2) gives no effect, rja = 1/2. The 



second term in the above equation represents Cooperon 
decay due to the probability of particle escape into the 
right hand reservoir. Multiple Andreev reflection at the 
boundary, illustrated in a sketch in Fig. 2b, causes a fac- 
tor of two reduction in the case of a polydomain ferro- 
magnetic wire connected to a superconducting reservoir. 
In the case of a ferromagnetic wire connected to an insu- 
lating reservoir, rja = 0, since the probability of particle 
escape into the right hand reservoir is totally supressed. 



In Fig. 4 we plot r]a for a ferromagnetic wire with arbi- 
trary polarization connected to a superconducting reser- 
voir. It is shown for each spin channel separately and we 
assume that the up-spins are the majority spins, (T| > (T[. 
For zero spin polarization, Jyt = ?7x = 1/2, as described 
above. On increasing the spin polarization, 77^ < 1/2 
whereas rj^ > 1/2. Majority carriers at the supercon- 
ducting interface cannot find minority carrier states to 
Andreev reflect into and are normally reflected which in- 
creases the probability to return whereas minority car- 
riers at the interface have an enhanced probability to 
escape. In the limit of total spin polarization, 77^ — )■ 
and rii^ —^ 1 which means that the reservoir appears to 
be totally insulating to majority spins but normal as far 
as minority spins are concerned. 
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FIG. 1. Schematic of the ferromagnetic wire comiected to a normal (N) reservoir on the left hand side, x 
and a superconducting (S) reservoir on the right hand side, x — 0. 
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FIG. 2. (a) Cooperon decay at the F/N boundary due to electron escape into the N reservoir, (b) Multiple Andreev 
reflection at the S reservoir changes the Cooperon boundary conditions at the interface. 
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FIG. 3. The prefactor, 7, of the increase, i5rf — Sr^, (Eq. (|l|)) in the weak locahzation contribution to the resistance 
as the reservoir changes from the normal state to superconducting, plotted as a function of the degree of polarization, 

for L/Ls — 20. The solid line is the total weak localization contribution whereas the dashed line does not include 
the contribution of density variations. The error increases as ^ — > 1 (see comments at the end of section III). 
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FIG. 4. The prefactor, rja, of the contact term in the spatiaUy averaged return probabiUty, 6Da, (Eq. (All)) for 
the F/S system, plotted as a function of the polarization, C, for L/Lg = 20. The solid line is for the majority spin 
channel whereas the dashed line is for minority spins. 
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